arXiv:1501.02060vl [physics.optics] 9 Jan 2015 
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We introduce a class of metamaterials with uniformly balanced gain and loss associated with com¬ 
plex permittivity and permeability constants. The refractive index of such a balanced pseudo-passive 
metamaterial is real. An unbounded uniform pseudo-passive medium has transport characteristics 
similar to those of its truly passive and lossless counterpart with the same real refractive index. 
However, bounded pseudo-passive samples show some unexpected scattering features which can be 
further emphasized by including such elements in a photonic structure. 


PACS numbers: 42.82.Et, 42.25.Bs, ll.30.Er 

In the last decade we have witnessed unprecedented 
advancements in realizing artificial materials which are 
specifically designed to exhibit features not found in na¬ 
ture. In the electromagnetic domain, such metamaterials 
are using their structural composition in order to obtain 
access of all four quadrants in the real e — pi plane. Var¬ 
ious exotic phenomena ranging from negative refraction 
to electromagnetic cloaking and super-lensing have been 
actively pursued w- While these opportunities have 
been extraordinary, optical materials exhibiting exotic 
values of permittivity 6 and/or permeability p are often 
prohibitively lossy. This is especially true for composite 
optical meta-materials. A natural solution to the prob¬ 
lem is to add a gain component and, thereby, to offset the 
losses. A pathway to achieve this goal has been recently 
proposed in Ref. |9j, and subsequently explored in a num¬ 
ber of publications QEH2S]. This proposal capitalizes on 
the notion of Parity-Time (VT) symmetry by utilizing 
balanced gain and loss elements which are judiciously 
distributed in space such that the complex index of re¬ 
fraction n(r) satisfies the relation n(r) = n*(—r). This 
research line gave many intriguing transport properties 
like double refraction, unidirectional invisibility, asym¬ 
metric transport, CPA/ Lasers etc. 

An alternative approach is to compensate the losses 
with gain while preserving the uniformity of the medium. 
This can be done simply by doping the lossy medium with 
active elements (dies), so that the doped material would 
have a real refractive index n(r) = y/epi = n*(r). The 
question, however, is whether such a balanced loss/gain 
medium with real and uniform refractive index will au¬ 
tomatically behave as a regular passive lossless medium. 
This is certainly true if the permeability pi is real, in 
which case the negative contribution to e" from gain will 
simply offset the positive contribution from absorption 

M. 

In this Letter we will concern with a different scenario 
where both 5 and pi are complex i.e. e = s' + is", pi = 
pi' + ipi". A typical example is plasmonic metamateri¬ 
als where the effective permeability is associated with 
the electric current induced in tiny ring resonators by 



FIG. 1: (color online) The electric E(z) (y-z plane) and the 
magnetic H(z) (x-z plane) field distribution of a resonance 
localized mode A r = 1.25 6pun for a Bragg grating with an em¬ 
bedded PP-defect (green). The grating consists of 20 lossless 
bilayers (orange and white) at each side of the defect. Their 
index of refraction is m = 1.45 and — 1.755 and their 
width is di = 0.2167 pan and gL = 0.179/ira. The pseudo¬ 
passive (PP) defect layer has index of refraction n p = 3.48 
corresponding to constituents e p — 3.48 — zO.Ol, pi p = e* while 
its width is d p = 0.0542/xra. 

oscillating magnetic fields. If both e and pi are com¬ 
plex, a real-value refractive index n can only be realized 
if £ = p 2 pi* (where p is real), which is quite possible in 
media with gain. In this case n = p\p\ and thus one 
expects that an unbounded medium consisting of such 
material will support traveling waves as its lossless pas¬ 
sive counterpart with the same index of refraction. We 
will therefore refer to media with real refractive index n 
and complex e" and pi" as pseudo-passive (PP) media. 
Surprisingly, once a pseudo-passive medium is confined 
in space, its scattering properties are completely different 
than the ones of its passive counterpart. The difference 
becomes even more striking when a pseudo-passive ma¬ 
terial is a component of a photonic structure. For exam¬ 
ple, the Fabry-Perrot resonances of a Bragg grating that 
contains pseudo-passive composite layers display sub- or 
super unitary transmittance at the band-edges. Remark¬ 
ably, under different circumstances, the same pseudo- 
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passive material can act as a perfect absorber, or it can 
trigger lasing instability. 

The scattering properties of a pseudo-passive struc¬ 
ture are understood better by considering monochro¬ 
matic wave propagation in an one-dimensional scattering 
set-up. In this case, the steady state (TE) electric field 
E(z ) is scalar and satisfies the Helmholtz equation 

0 (1) 

where co is the frequency of the field and c is the speed 
of light in vacuum while the spatially dependent index 
of refraction n(z) = e(z)ii(z ) is considered real for any 
position z. In the case that the pseudo-passive medium is 
extended over the whole space, Eq. 0 admits wave solu¬ 
tions of the form E{z) = E + exp (iknz) -\-E~ exp (—iknz) 
where k = uj/c is the free space wavevector. These solu¬ 
tions are identical to the ones found at an infinite lossless 
passive medium with index of refraction n. 

These similarities between a PP medium and a loss¬ 
less passive medium cease to exist once we turn to the 
analysis of the associated scattering problem involving 
bounded structures. For simplicity we assume that the 
domain which contains the pseudo-passive structure ex¬ 
tends in the interval 0 < z < L. We further assume that 
the scatterer is embedded in a homogeneous medium with 
uniform index of refraction no which, for simplicity, we 
consider it to be equal to unity i.e. no = 1. The solu¬ 
tion of Eq. 0 at the left of the scattering sample z < 0 
is El = E R exp (ikz) + E R exp (—ikz) while on its right 
z > L is Er = E^ exp(ikz) + E R exjp(—ikz). We can 
relate the amplitudes of forward and backward propa¬ 
gating waves on the left of the scattering domain with 
the amplitudes on its right via the transfer matrix M 


Er 

En 



( 2 ) 


In case of composite structures the transfer matrix (TM) 
M is a product of transfer matrices M n associated with 
each individual composite element i.e. A 4 = Yl n M n . 
The individual TM’s are evaluated by imposing the ap¬ 
propriate boundary conditions at the interface between 
consequent layers n and n + 1 of the structure. For a TE 
mode we have that the field itself is continuous while its 
spatial derivative satisfies the relation: 


"JL f dE n (z) \ 
Ji n \ dz J 


1 / dE n+1 (z) 

lhi +1 v dz 


interface 


(3) 


The transmission t and reflection r amplitudes for 
left and right incident waves are defined as El = 
Er/El ? r L = E L /Eft and tR = E L /E R , tr = E R /E R . 
They are written in terms of the TM elements as 


tL 


det M 
M-22 


, Lr = 


1 _ M-21 

M-22 ^ L M-22 


M 


12 


M 


22 


(4) 


Furthermore, we can show that det Ad = 1 and thus U, = 
tR = t. The associated transmittance T and reflectances 
Rl,r are T = \t\ 2 and 1Zl,r = \tl,r\ 2 - 

An alternative formulation of the scattering process 
is provided by the scattering matrix S which connects 
incoming and outgoing wave amplitudes i.e. 



The scattering matrix can be written in terms of the Ad- 
matrix elements as 


5 = 


r L tR \ _ 1 f —M2I 1 \ 

Tr J M22 V det(M) M\2 ) 


(6) 


and is useful for the theoretical analysis of lasing insta¬ 
bilities and perfect absorption. Below we provide some 
examples which illustrate the anomalous scattering prop¬ 
erties of a bounded structure that includes a PP-medium. 

We start our analysis with the investigation of the 
transport properties of a single PP-layer. The associated 
TM that describes the scattering process is 

M = D(L;no)~ 1 K(—; —)D(L;n p )K(—; —) (7) 
n 0 ii p n p /i 0 

where the matrices D and K are defined as 


k(LL,^ 

Tijn k*i 

D(L-n ) 


1 A + 

ni-(£)W i + (£)(W 

(exp {iknL) 0 \ 
y 0 exp {—iknL) J 


( 8 ) 


Using Eq. Q we evaluate the transmission and reflec¬ 
tion amplitudes of a single PP-layer: 


2 exp [—ikL\ 

2 cos[kn p L\—i(y p +^) sin [kn p L\ 

i (yp-^h) sin ( kn P L ) 

^ 2 cos [kn p L]—i(y p -\- sin [kn p L] 5 


(9) 


tr = t exp(— 2ikL) 


The subindex p indicates that we refer to the constituents 
of the PP-layer and y p = J jf- is its complex admittance. 


The Fabry-Perrot resonances are defined by the re¬ 
quirements T = 1 and R = 0 and occur at frequencies 
cjfp = kj?pc for which sm(kppn p L) = 0. This implies 
perfect resonant transmission with no losses and no gain, 
regardless of the value of y p and despite the fact that 
the constituents are complex. Let us now calculate the 
transport characteristics of the PP-layer at uj ^ ujpp. For 
example, at the middle between two consecutive Fabry- 
Perrot resonances, i.e. when cos [kn p L\ = 0, we have: 


r = 

A 


2 

2 

V — 

(Vp ^-/Vp) 

2/p + 1/ Up 

, /V — 
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_ 2 (yp Up) 
12/p + 11 2 


2 


= -8 
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( 10 ) 
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Remarkably, the off-resonance value of A is always neg¬ 
ative and independent of the slab thickness L indicating 
that its origin is associated with surface scattering rather 
than bulk scattering precesses. In fact, it is straightfor¬ 
ward to show that a PP-layer will result in amplification 
A < 0, for any frequency a;, irrespective of the sign of e". 

Next we consider the transport properties of a loss¬ 
less Bragg grating (BG) with one PP-defect layer. We 
assume, without loss of generality, that e" < 0. Before 
going on with our analysis we recall that a lossless pas¬ 
sive defect layer supports a resonant defect mode with a 
frequency uj r lying inside the photonic band gap of the 
BG. This resonance mode is localized in the vicinity of 
the defect layer and decays exponentially away from the 
defect. Near u; r , the entire composite structure displays 
a strong resonant transmission with T = 1 (and thus 
1Z == 0) due to the excitation of the localized mode. 


energy, due to the PP-layer, gives A ~ \Ed\ 2 lm(e") < 0 
and thus T + 7Z = 1 — A 1. The transmission and 
reflection spectrum are shown in Fig. and confirm the 
previous expectations. 

The opposite scenario is observed in the case that 
the magnetic field H(z) has a maximum in the domain 
around the defect layer. In this case T + 1Z = l — ic 1 
(we recall that A ~ \Hd\ 2 Tm{fi") > 0) and thus we have 
attenuation of the incident light. 

The amplification/attenuation mechanism is bet¬ 
ter investigated by introducing the overall amplifica- 
tion/absorption coefficient 0(e",u;) defined as the ratio 
of the total intensity of outgoing to incoming waves: 


e(e'» 


|£ZI 2 + |£AI 2 

i^i 2 + i^i 2 


(ii) 


The above expression can be written in the following form 



FIG. 2: (color online) (a) The transmittance T(oS), re¬ 

flectance 1Z(uj) and absorption A(uS) for the set-up of Fig. 
[l] (e 7/ < 0). A super-unitary transmittance and reflectance 
is evident in the frequency domain around the resonance fre¬ 
quency. (b) The amplification/absorption coefficient 0 versus 
frequency has a singular point at resonance frequency which 
signifies a lasing action, (c) The same set-up as in Fig. [I] for 
the case of CPA frequency. In this case the electric field as 
a minimum at the defect layer while the magentic field has 
a maximum (d) The amplification/absorption coefficient 0 
versus frequency shows a zero at the CPA frequency. 

A completely different scenario emerges when the de¬ 
fect layer consists of a PP-medium. We find that the res¬ 
onance localized mode has transmittance and reflectance 
values which are larger (smaller) that unity in distinct 
contrast to the case of lossless defect layers. The super 
(sub)-unitary transmission/reflection is related with the 
distribution of the electro-magnetic field around the PP- 
defect. Let us consider the case for which the electric field 
E{z ) takes its maximum value (anti-nodal point) at the 
position of the defect layer (see Fig. [I]). Obviously the as¬ 
sociated magnetic field H(z) will be having a nodal point 
at the defect layer. Then, an estimation of the absorbed 



where we have used the fact that det Ad = 1. The 
case 0(e",u;) > 1 indicates that an overall amplification 
has been achieved at the system. The opposite limit of 
0(e",u;) < 1 corresponds to attenuation. The two ex¬ 
treme cases of ©(e^,^) -A oc and 0 (cq PA ,cjcpa) 0 
indicate that the system has reached a lasing instability 
or behaves as a coherent perfect absorber (CPA) respec¬ 
tively. The condition for the former case is M. 22 ^" ^ w) = 
0 which is satisfied for some values e" = e'[ and uj = ujl- 
In fact, the complex zeros of Ad 22 correspond to the poles 
of the scattering matrix S', see Eq. lie" = 0 they lie 
at the lower part of the complex-cj plane due to causal¬ 
ity. As |e /7 | increases the poles move towards the real axis 
and at a critical value e" = one of these poles becomes 
real uj = ujr , thus signifying the transition to a lasing ac¬ 
tion. We have confirmed in Fig. Wp that the singularity 
in 0(e" = e A ,u; = uJl) —^ oc corresponds to the lasing 
point as it is calculated from the analysis of the poles of 
the associated scattering matrix. 

The other limiting case of CPA is achieved when the 
two terms in the numerator of Eq. ( [I2| ) become simul¬ 
taneously zero i.e. when E^ = E/kAd 2 i and E^ = 
—E^/M 12 . These two relations are simultaneously sat¬ 
isfied when M 21 M 12 T 1 = 0 —Y M 11 M 22 = 0 (re¬ 
call that det Ad = 1). The CPA frequency uj = cjcpa 
is evaluated from the condition Ad h(co’cpa) = 0 (the 
frequencies for which Ai 2 2 {^) = 0 are excluded since 
they correspond to lasing action). It should be stressed 
that a CPA requires coherent incident fields which sat¬ 
isfy an appropriate phase and amplitude relationship 
Er = Ad 21 = — -Ex/Ad 12 . An example of a CPA 

is shown in Fig. §i- 

The above super/sub-unitary scattering features be¬ 
come more pronounced when one of the composite layers 
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of the Bragg grating is substituted by a pseudo-passive 
medium. Specifically, the Fabry-Perrot (FP) resonances 
which are closer to the band edges will correspond to 
photons with small group velocities. Therefore each in¬ 
dividual photon at this frequency resides in the PP-layers 
for a long time. This long time interaction with the PP- 
medium leads to strong amplification/attenuation fea¬ 
tures. The magnitude of the enhancement/suppression 
is a growing function of the total thickness of the stack, 
as opposed to the case of a uniform pseudo -passive slab. 
In Fig. [3] we show a typical transmission spectrum of 
such composite structure. 



FIG. 3: (color online) The electric E(z) (blue line at y-z 
plane) and the magnetic H ( z ) (red line at the x-z plane) fields 
of two representative FP modes which reside at the upper (left 
inset) and lower (right inset) band edges. The grating consists 
of 50 bilayers (orange indicates a PP-layer and white a lossless 
passive layer). In the insets we show the field profiles for the 
first half of the structure, i.e. 25 layers. Left (right) inset 
correspond to the FP resonance at the upper (lower) edge of 
the lower (upper) band. The constituents of the PP-layers are 
e p = 4 — 0.001 z,AG — e p with a width layer d p = 0.1309 fim. 
The lossless passive layers have constituents e = 3 and g — l 
and width d — 0.3023 fim. The transmittance T(cj), of the 
grating is shown in the main panel. 

In the particular case of a FP resonance close to a pho¬ 
tonic band edge of a BG, consisting of two layers A and B, 
the field inside the periodic structure is a superposition of 
the forward and backward Bloch waves with Bloch wave 
number close to either the middle of the Brillouin zone, 
or to its boundary. In either case, the middle of each 
individual layer coincides with the node or the antinode 
of the respective standing wave inside this particular uni¬ 
form layer. The nodes (antinodes) of oscillating electric 
field coincide with the antinodes (nodes) of oscillating 
magnetic field and vice versa. 

At the same time, if the lower edge of a photonic band 
gap has electric field node (magnetic field antinode) in 
the middle of each A-layer, then the upper edge of the 
same photonic band gap has electric field antinode (mag¬ 


netic field node) in the middle of each A-layer [33] . As a 
consequence, Fabry-Perrot resonances close to neighbor¬ 
ing photonic band edges will have different dominating 
components - electric or magnetic - of the oscillating field 
inside a particular layer. 

In our example in Fig. [3] the A-layers consists of a PP- 
medium with negative e" and positive /x", which means 
that the Fabry-Perrot resonances with the dominant elec¬ 
tric field component in the A-layers will be enhanced, 
while the resonances with dominant magnetic field com¬ 
ponent in the A-layers will be suppressed. This is exactly 
what we see in Fig. [3J where the resonant transmission at 
every other photonic band edge is enhanced (suppressed). 

In conclusion, we have introduced a class of pseudo¬ 
passive metamaterials for which the losses are balanced 
by gain in a way that the index of refraction is real and 
uniform throughout the medium. Although the light 
propagation in an unbounded pseudo-passive medium 
has the same characteristics as in a passive lossless 
medium with the same index of refraction, their scat¬ 
tering properties differ dramatically. When such me¬ 
dia are incorporated in a photonic structure, it can lead 
to super/sub-unity transmittance/reflectance indicating 
strong absorption and amplification mechanisms for the 
total structure. In many occasions these effects can turn 
to a coherent perfect absorption of incident waves or to 
lasing instabilities. It will be interesting to investigate 
the realization of these phenomena under vectorial con¬ 
ditions and in higher dimensions. 
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